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II 

REVIEW  OF  RESEARCH  COMPLETED 

The  general  objective  of  the  research  effort  has  been  a  better  definition 
of  the  explosion  and  earthquake  source  processes.  Specific  elements  of  the 
research  program  are:  1)  recording  of  broadband  data  from  events  at  the 
Nevada  Test  Site;  2)  analysis  of  the  coherence  of  ground  motion  near  explosions 
and  earthquakes;  3)  study  of  the  relative  isotropic  and  non-isotropic  components 
of  explosive  sources  through  the  application  of  moment  tensor  inversion 
techniques;  4)  analysis  of  regional  surface  wave  data  in  order  to  obtain 
models  for  the  velocity  and  attenuation  of  the  crust;  and  5)  archival  of 
near  and  regional  data  sets  which  are  of  value  to  the  general  discrimination 
problem. 

Table  1  contains  a  list  of  research  contributions  which  have  been  completed 
during  the  grant  period.  All  of  these  have  been  published  except  for  items 
3  and  4.  Item  3  is  a  condensed  version  of  material  that  appeared  in  our 
Technical  Report  for  01  October  1982  -  31  October  1983,  so  it  is  not  reporduced 
here.  Item  4  is  included  as  Section  III  of  this  report. 

During  the  grant  period  two  explosions  at  NTS  were  recorded  with  our 
digital  network.  The  Coalora  event  was  detonated  in  Yucca  Flat  on  February 
11,  1983.  It  was  recorded  by  10  three-component  accelerometer  stations  in 
the  distance  range  0.7-5. 4  km  and  an  azimuthal  range  of  150  degrees.  The 
Chancellor  event  was  detonated  in  Pahute  Mesa  on  September  1,  1983.  It 
was  recorded  by  11  three-component  accelerometer  stations  in  the  distance 
range  1.8-9. 7  km  and  an  azimuthal  range  of  360  degrees. 
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THE  TIBETAN  LITHOSPHERE: 

NEW  SEISMIC  DATA  FAVOR  AN  OLD  MODEL 


A in-  Yip  Chun 
Thomas  V.  McEuilly 

Seismographic  Station 
and 

Earth  Sciences  Division 
Lawrence  Berkeley  Laboratory 

University  of  California 
Berkeley,  California  94720 


Introduction 

The  continuing  convergence  of  the  Indian  and  Eurasian  continents 
since  the  Eocene-Oligocene  transition  some  40  Myr  ago  has  produced  a 
zone  of  deformation  extending  as  much  as  3000  km  northeast  of  the 
Himalayan  range.1  Rising  to  5  km  average  elevation  over  an  area  almost 
the  size  of  the  entire  U.S.  western  cordillera,  the  Tibetan  Plateau  is  one 
of  the  most  conspicuous  land  features  on  this  planet.  Geophysicists 
debate  whether  its  underlying  70-km-thick  crust2-6  was  created  by 
underthrusting  of  the  Indian  subcontinent7,8  (the  underthrusting  model), 
by  horizontal  shortening  with  thickening  under  N-S  compression9  (the 
contraction  model),  or  by  some  other  process. 

Chinese  and  French  scientists  recently  reported10  new  paleomag- 
netic  data  suggesting  that  the  Tibetan  landmass  is  made  up  of  several 
continental  and/or  island  arc  fragments  which  were  progressively 
attached  to  the  Eurasian  continent  between  late  Paleozoic  and  late  Cre- 
taceous11,12  (the  accretion  model).  Crustal  shortening  by  the  associated 
underthrusting  brings  the  crustal  thickness  to  70  km,  twice  the  normal 
continental  value.  The  contraction  and  accretion  models  emphasize 


are  known. 


Tibetan  crustal  velocity  models  have  been  derived  mainly  from 
seismic  Rayleigh  and  Love  surface  wave  dispersion  data.2-6  Recent 
investigations17,18  have  quantified  model  parameter  resolutions  with  for¬ 
mal  geophysical  inversion  techniques  concluding,  generally,  that  the 
Tibetan  crust  is  about  70  km  thick,  with  a  V5  in  the  lower  half  crust 
averaging  3.60  km/sec,  a  value  sufficiently  low  to  be  considered  indica¬ 
tive  of  partial  melting  5  These  surface  wave  studies  all  use  mixed-path 
data,  i.e.,  wave  propagation  paths  sampling  both  Tibetan  and  neighboring 
crust,  in  which  the  fraction  of  Tibetan  path  rarely  exceeds  2/3  the  total 
path.  Numerical  simulations  readily  demonstrate  that  formal  resolution 
analyses  fail  when  plane-layered  models  are  applied  to  a  mixed- 
structure  having  large  lateral  variations  in  properties. 

Refraction  studies5  yield  high  upper  mantle  VP  and  Vs  values  of 
8.12  ±  0.06  km/sec  and  4.80  ±  1.0  km/sec,  respectively.  High  frequency 
Sn  (0.5  to  2  Hz),  a  seismic  phase  considered  diagnostic  of  physical  pro¬ 
perties  in  the  uppermost  mantle,  is  found  to  propagate  efficiently 
beneath  Tibet  and  the  Indian  Shield.19  This  is  in  marked  contrast  to  poor 
Sn  propagation  seen  in  regions  of  high  heat  flow  and  recent  volcanism, 
such  as  the  western  U.S.  Basin  and  Range  province,  the  East  African  Rift 
system20  and  the  northern  Iranian  Plateau.21 

Taken  together,  these  results  from  previous  surface  wave  and  body 
wave  studies  imply  an  intriguing  phenomenon  in  which  a  thickened,  hot 
Tibetan  crust  apparently  overlies  a  cold,  shield-like  mantle. 
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invariant  from  one  propagation  path  to  the  other.  Slowness  data  in  the 
8-100  sec  period-range  were  fitted  by  least-squares,  as  shown  in  Fig.  3a, 
for  each  period.  The  reciprocal  of  the  slowness  intercept  gives  the 
desired  pure-path  Tibetan  group  velocity.  Figures  3b  and  3c  summarize 
the  resulting  pure-path  and  several  representative  mixed-path  group 
velocity  curves  for  Rayleigh  and  Love  waves.  A  TPF  of  0.5  (50%  Tibetan 
path)  is  typical  for  the  data  used  in  most  previous  surface  wave  studies 
of  Tibet.  From  Figures  3b  and  3c,  relative  to  our  estimated  pure-path 
velocities,  the  non  regionalized  group  velocities  are  in  error  by  about  0.4 
and  0.2  km/sec,  respectively,  for  Rayleigh  and  Love  waves.  Inversion  of 
such  mixed-path  data  will  lead  to  serious  errors  in  the  resulting  struc¬ 
ture. 

The  validity  of  our  regionalization  scheme  can  be  tested  with 
independent  regionally  recorded  mixed-path  data  sets,  by  the  ability  of 
the  standard  curves  in  Figures  3b, c  and  the  Tibetan  boundary  in  Figure  1 
to  predict  the  independent  data  set.  To  illustrate,  we  consider  the  con¬ 
veniently  tabulated  mixed-path  group  velocity  data  for  Rayleigh  and  Love 
waves  from  a  previous  study4  involving  17  paths  along  which  the  average 
TPF  is  0.66,  possibly  the  largest  of  any  previous  Tibetan  studies.  When 
the  results  in  Figures  3b, c  are  used  to  construct  Rayleigh  and  Love  wave 
standard  curves  with  a  TPF  of  0.66,  they  are  found  to  fit  the  correspond¬ 
ing  mixed-path  data4  quite  closely,  with  a  root-mean-squares  difference 
of  about  0.03  km/sec.  The  velocity  model  determined  from  the  TPF=0.66 
mixed-path  data  is  shown  in  Figure  4  labeled  TP4. 
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In  the  bottom  40  km  of  the  Tibetan  crust,  where  previous  mixed-path 
studies  have  persistently  reported  an  anomalously  low  3.6  km/sec  aver¬ 
age  velocity,  we  find  Vs  of  3.88  ±  0.08  km/sec,  which  compares  with  the 
3.85  km/sec  value  found  in  the  lower  crust  of  the  Canadian  Shield.32 

Our  surface  wave  data  do  not  constrain  details  in  the  upper  mantle 
beneath  the  thick  Tibetan  crust;  the  inversion  yields  a  4.50-4.55  km/sec 
average  upper  mantle  Vs,  a  value  in  qualitative  agreement  with  analyses 
of  travel-time  residuals  which  indicates  that  the  mean  velocity  in  the 
topmost  200-300  km  of  the  Tibetan  mantle  must  be  anomalously  low.5,33 
This  contrasts  with  the  high  8.12  km/sec  Vp  and  4.80  km/sec  Vy  obtained 
under  Tibet  in  mantle  refraction  studies.5  Taken  together,  these  results 
imply  a  thin,  high-velocity  mantle  lid  extending  no  deeper  than  about 
100  km  beneath  the  surface.5 

Our  proposed  crust-mantle  shear  velocity  model  TP  100  of  the 
Tibetan  Plateau  is  shown  in  Figure  4,  along  with  the  calculated  and  <he 
observed  surface  wave  data.  By  normal  standards  the  fit  is  excellent. 
Velocity  uncertainties  are  given  for  the  layer  thicknesses  shown. 


Temperature  at  Base  of  Crust 

Recent  research34  yields  a  correlation  between  the  upper  mantle 
compressional  velocity  Pn  and  estimated  temperature  Tm  at  the  base  of 
the  crust.  The  equation  Pm  =  8.456  -  0.000729  Tm,  where  Pm  is  the  meas¬ 
ured  Pn  corrected  to  a  35  km  reference  depth  using  a  pressure  deriva¬ 
tive  of  0.015  km/sec  per  100  MPa,  is  based  on  an  exhaustive  study  cover¬ 
ing  numerous  distinct  tectonic  provinces  across  North  America.34  Using 
8.12  ±  0.06  km/sec  as  the  Tibetan  Pn,  this  equation  gives  a  Tm,  at  the 
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develop  under  shear  in  the  lowest  part  of  the  overlying  Eurasian  crust, 
and  water  forced  from  the  uppermost  part  of  the  underthrusted  Indian 
continental  crust  enhances  the  cracks  by  hydraulic  fracturing  or  by 
chemical  corrosion.38  Pore  pressure  then  rises  as  a  direct  result  of 
shear-straining  and  the  local  heating  resulting  from  it. 

Beneath  the  Himalayas  and  southern  Tibet,  where  the  lower  crust  is 
recently  emplaced,  the  amount  of  water  being  released  should  be  rela¬ 
tively  significant,  lowering  the  sliding  friction  in  the  fault  zone.  Some  of 
the  water  may  find  its  way  to  the  surface  through  fissures,  creating  hot- 
springs  found  in  that  region.  Beneath  northern  Tibet,  in  contrast,  lies 
the  dry,  long-heated  leading  edge  of  the  underthrusted  Indian  litho¬ 
sphere,  where  a  combination  of  pre-heating  and  increased  sliding  fric¬ 
tion  may  cause  local  partial  melting.  This  south-to-north  gradation  can 
explain  both  the  only  active  volcanism  (Fig.  1)  and  the  near  absence  of 
geothermal  activity  in  the  Kun-Lun. 

A  Missing  Granitic  Layer? 

The  uniformly  high,  3.88  km/sec  V5  in  the  40-km-thick  lower  crustal 
section  does  not  appear  compatible  with  the  continental  underthrusting 
model.  Apparently  lacking  is  a  typical  granitic  layer  (V5  ~  3.65  km/sec) 
expected  in  its  upper  portion.  A  case  can  be  made,  however,  for  its 

upper  half  being  a  granitic  layer  in  disguise. 

»  1  % 

As  a  result  of  underthrusting,  the  average  lithostatic  pressure  in  the 
granitic  layer  (assuming  a  20  km  thickness)  of  the  Indian  crust  is 
increased  from  300  to  1400  MPa;  and  the  temperature  to  no  more  than 
500°  C  in  the  available  time  since  large-scale  underthrusting  commenced 
in  Miocene  time.  Mean  pressure  and  temperature  derivatives  of  0.0194 
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discriminate  among  competing  models  for  the  development  of  the 
Tibetan  Plateau.  Major  results  are: 

1.  The  crust  beneath  Tibet  is  74  ±  10  km  thick  with  average  Vs  of  3.54 
km/sec. 

2.  A  prominent  shear  wave  low  velocity  zone  is  seen  in  the  depth  range 
of  24-34  km,  where  V5  is  approximately  2.64  km/sec. 

3.  V5  in  the  lower  40  km  of  the  crust  is  uniformly  high  (3.88  ±  0.08 
km/sec). 

4.  The  LVZ,  overlain  by  an  upper  crust  with  properties  similar  to  the 
Basin  and  Range  or  the  Andes,  and  underlain  by  a  high-velocity, 
presumably  cool,  lower  crust,  is  better  explained  by  high  pore  water 
pressure  than  by  partial  melt. 

5.  The  high  V5  values  in  the  lower  crust,  along  with  previous  body  wave 
studies  of  upper  mantle  Pn  and  5^,  imply  the  presence  of  an  unusu¬ 
ally  cool  region  in  the  depth  range  50-100  km. 

These  results  may  be  interpreted  to  indicate  that  the  Eurasian 
geotherm  was  probably  similar  to  that  of  the  Basin  and  Range  until 
large-scale  continental  underthrusting  began  during  early  Miocene  time. 
The  underthrusting  Indian  lithosphere  some  120  km  thick43  was  warmed 
from  above  by  the  heat  generated  in  the  fault  zone  and  from  below  by  the 
heat  flux  from  the  underlying  asthenosphere.  The  short  time  available 
does  not  allow  thermal  equilibration  near  its  mid-section. 

We  present  the  schematic  drawing  in  Fig.  5  to  describe  our  tectonic 
model  of  Tibet  consistent  with  the  above  conclusions  on  the  evolution  of 
the  plateau. 
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Figure  Captions 

Figure  1.  Surface  wave  paths  used  to  study  the  crust-mantle  structure 
beneath  the  Tibetan  Plateau  (shaded  region).  The  epicenters 
of  earthquakes  used  in  this  study  are  indicated  as  small  filled 
circles,  and  WWSSN  stations  are  shown  as  open  circles  with 
code  letters.  Active  volcanoes  are  marked  as  larger  filled  cir¬ 
cles.  Map  based  on  Molnar  and  Tapponnier.1 

Figure  2.  Long-period  surface  wave  seismograms  (above)  and  the 
results  of  the  moving-window  group  velocity  analysis  (below) 
for  a  0.75  fraction  of  Tibetan  Plateau  path  (TPF)  and  a  0.27  TPF 
path,  left  (Rayleigh  Wave)  and  right  (Love  Wave),  respectively. 
Limits  of  the  vertical  bars  and  shading  represent  two  different 
energy  contour  levels. 

Figure  3.  Pure-path  group  velocity  determination.  Figure  3a  demon¬ 
strates  for  two  periods  the  least-squares  determination  of 
pure-path  slowness  values,  the  intercepts  which  are  shown 
with  their  uncertainties  (two  standard  deviations).  Slownesses 
are  plotted  for  a  given  period  as  a  function  of  non-Tibetan  path 
percentage.  Figures  3b  and  3c  show  the  resulting  pure-path 
Rayleigh  and  Love  wave  group  velocity  curves  so  determined, 
along  with  some  mixed-path  curves. 
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Introduction 

VELA-sponsored  research  at  Berkeley  began  around  1960  with  the 
program  of  Dr.  Don  Tocher  and  Prof.  Perry  Byerly  to  expand  their  central 
California  network  of  the  UC  Seismographic  Stations,  and  to  link  for 
the  first  time  such  a  network  by  FM  telemetry  to  a  central  recording 
site  (see  Figure  1).  Their  purpose  was  to  improve  research  capabilities 
for  studying  the  mechanisms  of  local  earthquakes.  Shortly  thereafter, 
with  AFOSR  support,  the  BKS  WWSSN  station  was  installed,  and  a  three- 
component,  broadband  (flat  velocity  response,  0.03-10  Hz,  dual  gains) 
station  was  installed  at  BRK  with  continuous  slow-speed  (0.06  ips) 

FM  recording  on  magnetic  tape.  In  1964  the  tape  recording  system  was 
expanded  to  handle  also  six  of  the  short-period  telemetered  stations. 

During  the  subsequent  two  decades  leading  to  today,  the  AFOSR- 
supported  research  program  has  concentrated  on  investigations  into  the 
source  mechanisms  of  earthquakes  and  explosions,  along  with  techniques 
for  discriminating  between  them.  From  the  start,  the  research  was 
based  largely  upon  the  new  broadband  data  being  acquired  by  the  network 
from  local  earthquakes  and  from  underground  explosions  at  the  nearby 
Nevada  Test  Site  (NTS),  although  some  studies  used  teleseismic  data. 
Promising  results  of  early  discrimination  studies  with  NTS  data,  along 
with  source  studies  of  moderate-size  central  California  earthquakes, 
led  to  further  expansion  of  the  broadband  stations  to  BKS,  JAS,  WDC 
and  SAO  (see  Figure  2).  Early  work  on  near-field  observations  of 
earthquakes,  at  Berkeley  and  elsewhere,  led  to  initiation  in  the  early 
seventies  of  the  AFOSR  "Near  Field  Program",  a  multi-institution  co¬ 
operative  study  of  the  Bear  Valley  -  Stone  Canyon  seismogenic  zone 


Figure  1.  UC  Seismographic  Stations  network  in  1S64  with  third 

quarter  seismicity.  Stations  of  AFOSR-supported  telemetered 
network  are  underlined. 
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Spectral  Evidence  for  Fault  Rupture  Parameters 

The  5.5  1966  Parkfield  earthquake  was  well  recorded  on  the  BRK 
broadband  system,  providing  data  for  many  investigation  since.  In  one 
study  the  Love  wave  spectra  were  analyzed  for  propagation  of  the 
rupturing  fault  surface.  This  early  application  of  directivity  effects 
compared  Love  wave  spectra  for  five  of  the  larger  earthquakes  in  the 
sequence,  shown  in  Figure  1.  The  spectra,  seen  in  Figure  2,  reveal  clear 
notches  at  periods  consistent  with  a  fault  rupture  length  of  about  30  km 
and  a  rupture  velocity  of  2.2  km/sec.  The  four  smaller  events  (Ml  3.8 
to  5.1)  are  quite  similar,  scaling  at  10-20  sec  periods  closely  with  Ml- 
The  5.5  Ml  raainshock  at  these  periods  scales  to  a  much  larger  value 
than  its  measured  ML,  as  seen  in  Figure  2,  which  represents  an  early 
demonstration  of  the  ML  saturation  phenomenon. 
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Short  Period  Discriminant  at  Regional  Distances 

A  spectral  ratio  method  was  shown  to  successfully  discriminate 
between  underground  explosions  and  natural  earthquakes  when  applied  to 
the  Pg  phase  at  regional  distances.  The  data  set  included  69  events 
recorded  at  the  Berkeley  network  station  JAS  in  the  distance  range 
250  to  500  km  and  the  magnitude  range  2.8  to  4.5  (Table  1).  The 
Pg  spectral  ratio  (0.6-1.25  Hz)/ (1.35-2.0  Hz)  shows  explosions  and 
natural  earthquakes  to  separate  into  distinct  sets  for  magnitude 
greater  than  3.2,  with  explosions  relatively  richer  in  the  high  frequency 
band.  An  interesting  result  was  that  the  spectra  of  afterevents  of 
large  explosions  resemble  the  spectra  of  explosions  more  than  that  of 
natural  earthquakes.  However,  these  afterevents  appear  to  be  more  like 
natural  earthquakes  when  the  m^:!^  discriminant  is  applied. 
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Table  1  ( continued ) 


No. 

Type 

Date 

OT  (GMT) 

Lat  (N) 

Long  (W) 

A  (km) 
from  JAS 

M14> 

Spectral 
ratio  (/*,) 

Comments 

54 

c 

1968  December  12 

15 

20 

00 

<37  1) 

(1160) 

(390) 

3  8 

I  36 

(4) 

55 

c 

1968  December  19 

17 

30 

22-8 

37-2 

116-5 

360 

3-4 

1-24 

<:i — BENHAM  'Aftcrevent' 

56 

c 

1968  December  19 

19 

18 

19-6 

37-3 

1 16-4 

360 

3  5 

0-48 

"'—BENHAM  'Afierevent' 

57 

c 

1968  December  19 

19 

54 

01-2 

37  2 

116  5 

360 

3  4 

0-39 

<J’— i BENHAM  •Aftercvcnf 

58 

c 

1968  December  19 

22 

23 

26-3 

37-2 

116  5 

360 

3-6 

1  06 

"’ — BENHAM  'Afierevent'  (collapse? 

59 

c 

1968  December  20 

20 

08 

20-4 

37-2 

116  5 

360 

3  8 

0  61 

"’ — BENHAM  'Afierevent' 

60 

c 

1968  December  21 

00 

14 

25-2 

37-3 

1 16-5 

360 

4-3 

— 

*” — BENHAM  'Afterevent) 

61 

c 

1969  January  6 

06 

34 

14-5 

37  3 

1 16-5 

360 

4-2 

— 

"’ — BENHAM  'Aftcrevent' 

62 

c 

1969  January  10 

09 

41 

21-5  * 

37-2 

116  5 

360 

3-9 

1-72 

"' — BENHAM  'Afterevent' 

63 

c 

1969  January  10 

17 

01 

44-5 

37-2 

116  5 

360 

3-8 

2-82 

"’—BENHAM  ' Afterevent' 

64 

c 

1969  January  10 

17 

14 

17  2 

37-2 

116  5 

360 

3  8 

0  88 

"’ — BENHAM  ‘Afterevent' 

65 

c 

1969  March  18 

14 

40 

02-7 

37-2 

116-0 

400 

3-8 

069 

in 

66 

c 

1969  September  16 

IS 

43 

49-2 

37-2 

116  5 

360 

3-8 

2-38 

— JORUM  'Afierevent' 

67 

c 

1969  September  16 

16 

23 

53-8 

37-3 

116  5 

360 

3  9 

1-39 

"'—JORUM  'Aftcrevent' 

68 

c 

1969  September  16 

17 

31 

14-7 

37-3 

116-5 

360 

40 

3  44 

"'—JORUM  'Afierevent' 

69 

c 

1969  September  16 

18 

15 

39-3 

37-3 

116  5 

360 

3  8 

1  52 

— JORUM  ‘Afterevent' 

111  Magnitudes  based  on  P,  and  P,  amplitudes  at  JAS.  The  scale  is  based  on  Wood-Anderson  magnitudes  at  BRK  for  larger  events. 

Personal  communication.  Don  Springer,  1970.  The  letter  designations  (Y)  and  (P)  indicate  the  test  areas  Yucca  Flat  and  Pahute  Mesa,  respective!) ; 
(R)  indicates  an  elongated  explosive  charge  was  detonated  as  a  ‘  row  shot  ’  at  Buckboard  Mesa.  The  origin  times  of  collapses  may  be  uncertain  b>  a  feu 
seconds. 

U.S.C.G.S.  hypocentre. 

Arrival  times  at  the  Berkele>  net  indicate  NTS  hypocentre.  A  =  360  or  390  km  is  assumed  for  purposes  of  attentuation  and  magnitude  calculations 

*  Arrival  times  available  for  this  event  are  inconsistent  so  that  the  epicentral  distance  is  uncertain.  The  magnitude  thus  may  be  as  large  as  3-2  if  the  epicentre 
is  500  km  from  JAS.  The  error  in  the  attenuation  correction  due  to  the  epicentral  uncertainty  is  not  significant. 


*  -'NcM'ol'  Ear "iquoV.cs 

*  -  F (Close's 

*  -  Explosion  ’eftcreveets’ 


f.'cj-vTviae 


Fig.  I  P,  spectral  ratio  of  vertical  surface  displacement,  corrected  for  attenuation 
with  Q  —  400,  at  JAS  for  events  within  100  km  of  NTS.  I  and  II  denote  the 
earthquake  and  explosion  fields,  respectively. 
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Seismograpnic  Stations 


Station 

Type  of  instrument 

Belleview,  Florida 
(BEFL) 

Geotech  18300  (CANNIKIN) 

Berkeley,  California 
(BRK) 

Press  Ewing  LP  (CANNIKIN) 

Houlton,  Maine 
(HNME) 

Portable  Benioff  (MILROW) 
Geotech  18300  (CANNIKIN) 

Kanab,  Utah 
(KNUT) 

Large  Benioff  (MILROW) 
Geotech  18300  (CANNIKIN) 

Las  Cruces,  New'  Mexico 
(LCNM) 

Portable  Benioff  (MILROW) 

Priest,  California 
(PRI) 

Portable  Benioff  (CANNIKIN) 

Red  Lake,  Ontario 
(RKON) 

Portable  Benioff  (MILROW) 
Geotech  18300  (CANNIKIN) 

San  Jose,  Texas 
(SJTX) 

Geotech  18300  (MILROW) 
Geotech  18300  (CANNIKIN) 

Location 

Epicentra!  dist. 
from  MILROW 
(deg) 

Azimuth 
from  MILROW 
(deg) 

28°  54'  19' N; 

82  03  52  W 

73  1 

64-8 

37  52  -4N; 

122  15  -6W 

42-7 

84-8 

46  09  43  N; 

67  59  09  W 

66-9 

44  0 

37  01  22  N; 
112  49  39  W  - 

49- 1 

79-1 

32  24  08  N; 
106  35  58  W 

56- 0 

79- 1 

36  08  -5N; 

120  39  -9W 

44-8 

85-6 

50  50  20  N; 

93  40  20  W 

51*5 

540 

27  36  43  N; 

98  18  46  W 

64-4 

77-2 

Sec onds 


Fig.  I  Mean  impulse  train  (solid  line)  ±  its  standard  deviation  (dashed  line) 
for  MILROW.  Impulse  trains  from  KNUT,  HNME,  LCNM  and  RKON  were 
used  to  compute  the  mean. 


.>? :  i-i', 


Fig.  2  Mean  impulse  train  (solid  line)  ±  its  standard  deviation  (dashed  line)  for 
CANNIKIN.  Impulse  trains  from  HNME,  KNUT,  SJTX  and  BRK  were  used 
to  compute  the  mean. 
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Unidentified  NTS  events  listed  with  explosions  on  basis  of  Rayleigh-/1,,  ratios. 
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:  laboratory  seismic  network 


-  network  (open  triangles)  and  the  Berkeley  (BRK) 


Broadband  Studies  of  Earthquakes  in  the  Near  Field 

Broadband  recordings  (0.03  to  10  Hz)  in  the  near  field  (2  to  40  km) 
were  used  to  study  a  series  of  earthquakes  along  the  San  Andreas  fault 
in  central  California.  Special  broadband  instruments  were  operated  on 
both  sides  of  the  fault  at  the  San  Andreas  Geophysical  Observatory  (SAGO). 
The  13  earthquakes  used  in  the  study  are  listed  in  Table  1  and  plotted 
in  Figure  1.  Spectra  were  used  to  estimate  scalar  moments  (Figure  2) 
and  establish  a  relationship  between  moment  and  local  magnitude  (Figure  3). 
The  near  field  terms  of  the  elastic  displacement  field  caused  by  an 
earthquake  were  clearly  present  on  the  seismograms  (Figure  5)  and  at 
near  stations  could  be  reasonably  well  modeled  by  a  point  dislocation 
source  in  a  homogeneous  halfspace  (Figures  4  and  5),  although  the  effects 
of  tilt  on  the  seismometers  also  had  to  be  included. 
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Fig.  2  Example  illustrating  two  different  estimates  of  the  low-frequenc>  level  fi0  and  the  corner 
frequents  /„  of  the  spectrum  from  the  NS  component  at  SAGO-East  from  esent  6  (Table  I). 


Fig.  3  Relation  between  seismic  moment  and  magnitude.  The  solid  line  was  fit  to  the  data  points  h> 
linear  regression  and  the  dashed  line  is  from  W>ss  and  Brune  (l%8t  The  SAGO-East  values  have  been 
divided  by  a  factor  of  3  to  remove  an  amplification  effect. 
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Near  Field  Project 

The  "Near  Field  Project"  grew  out  of  some  of  t.he  problems  that  became 
apparent  at  the  Woods  Hole  meeting  on  seismic  discrimination  in  1970. 
Because  the  type  of  high  quality  data  necessary  to  test  various  theoretical 
models  of  an  earthquake  was  not  available  at  that  time,  a  cooperative 
experiment  was  designed  to  trap  a  moderate  size  earthquake  within  a  network 
of  stations  designed  to  provide  data  over  a  large  range  of  frequency, 
azimuth,  and  distance.  The  Stone  Canyon-Bear  Valley  section  of  the  San 
Andreas  fault  in  central  California  was  selected  as  the  target  area, 
and  U.C.  Berkeley  had  the  responsibility  of  developing  and  installing 
a  network  of  three-component  broadband  seismographs.  Both  acceleration  and 
displacement  were  recorded  and  an  effective  bandwidth  of  0.02  to  50  Hz 
was  achieved  (Figure  2). 

The  9  stations  of  the  network  were  installed  in  early  1973  and 
remained  operational  until  early  1977.  Figure  1  shows  the  locations  of 
the  stations  and  the  epicenters  of  3  earthquakes  which  provided  useable 
data.  Figure  3  shows  the  recordings  of  ground  acceleration  from  the 
magnitude  3.3  event  on  July  6,  1974,  obtained  at  an  epicentral  distance 
of  2  km.  The  clear  separation  of  P  and  S  waves,  the  short  duration  of 
the  phases,  the  high  frequencies,  and  the  rather  large  accelerations 
for  an  event  of  this  size  are  all  of  interest.  Figure  4  shows  the 
recordings  of  ground  displacement  as  the  same  station  from  the  same 
event.  Here  we  see  the  large  effect  which  can  be  introduced  by  ground 
tilt  at  stations  in  the  near  field. 


FULL  SCALE  INPUT  ACCELERATION 
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Near  Field  Experiments  at  the  Nevada  Test  Site 

Beginning  in  1969,  U.C.  Berkeley  has  conducted  an  experimental  program 
of  recording  explosions  at  the  Nevada  Test  Site.  Table  1  lists  the  13 
events  which  have  been  recorded  so  far.  These  experiments  have  included 
events  in  both  Yucca  Valley  and  Pahute  Mesa.  They  have  also  included 
a  variety  of  recording  arrangements,  ranging  from  arcs  at  a  single  distance, 
networks  containing  stations  at  a  variety  of  azimuths  and  distances, 
one-dimensional  arrays,  and  two-dimensional  arrays.  Figure  1  shows  the 
events  and  stations  for  the  experiments  at  Pahute  Mesa.  The  combined 
data  set  from  this  series  of  experiments  is  now  large  enough  to  permit 
systematic  studies  of  propagation  and  site  effects,  particularly  for  the 
Pahute  Mesa  region. 
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NTS  Area  Crustal  Velocity  and  Q  Structure  from  Broadband  Surface  Waves 
An  application  of  the  phase-matched  filtering  process  to  multimode 
Rayleigh  and  Love  waves  from  NTS  explosions  recorded  at  the  four  LLNL 
broadband  stations  (Figure  1,  Table  1)  proved  surprisingly  effective  in 
extracting  pure-path  group  velocities.  Fundamental  mode  Love  and  Rayleigh 
waves  and  first  shear  mode  of  the  Rayleigh  wave  were  observed,  yielding 
data  in  the  period  range  2  to  50  seconds.  Data  for  the  Landers  station 
are  shown  in  Figure  2,  in  which  observations  are  given  as  dots.  These 
velocities  are  the  lowest  found  for  the  Basin  -  Range,  apparently  because 
they  represent  the  shortest  observation  paths  ever  used,  thus  avoiding 
upward  bias  of  true  velocities  inherent  in  the  averaging  process  which 
is  involved  in  using  long  paths.  Inversion  for  structure  yields  the 
NTS-Landers  model  of  Figure  3,  and  the  calculated  dispersion  curves  are 
shown  as  solid  curves  in  Figure  2.  Statistically  significant  differences 
in  models  for  the  four  paths  are  seen,  apparently  correlating  with  the 
regional  variations  in  heat  flow.  Pure-path  amplitudes  are  also  extracted 
by  the  phase-matched  filtering,  and  they  allow  estimation  of  a  regional  Q 
structure  which,  while  highly  uncertain,  indicates  Q  values  greater  than 
100  in  the  upper  5-10  km,  and  less  than  50  below  that  depth.  These  new  dat 
on  Basin  -  Range  structure  will  be  used  to  improve  Greens  Functions  for 
explosion  source  inversion  from  regional  data  sets. 


SHEAR  VELOCITY  IN  KM/SEC 


Observed  (dots)  and  calculated  dispersion.  Figure  3.  NTS  -  Landers  model  based  upon  invers 

of  data  in  Figure  2. 
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Figure  3. 

Estimate  of  the  isotropic  part  of  the  first-degree  moment  tensor  of  the 
explosion  HARZER. 
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